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Abstract 



S^ ■ In this paper we study the interaction between the general form of viscous varying 

[~^ . modified cosmic Chaplygin gas and the Tachyon fluid in the framework of Einstein 

00 \ gravity. We want to reconstruct the Tachyon potential and total equation of state 

parameter graphically by using numerical methods. In the presence of deceleration 
■ ■ parameter, the interaction between components becomes sign changeable to explain 

(^ \ different evolutionary eras in the universe. We review the potential and total equation 

cn ' of state parameter in Emergent, Intermediate and Logamediate scenarios of scale fac- 

tor numerically. Analysis of total equation of state parameter show that, oJtot < — 1 
and ujtot > — 1 imply the phantom-like and quintessence-like behaviors respectively, we 
^ I have checked the effects of cosmic and viscosity elements on the interaction process. 

H ■ Stability is checked in all the models by the squared velocity of sound. 
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1 Introduction 

The most attractive subject in cosmology is the accelerating expansion of the universe which 
is based on the recent astrophysical data explaining the universe is spatially flat and an 
invisible cosmic fluid called dark energy with a hugely negative pressure which is responsible 
for this expansion. This type of matter violates the strong energy condition, i.e.,p + 3p < 0. 
There are various phenomenological models describing dark energy. Cosmological constant 
is the simplest one which gave rise to the A-CDM model, but it suffers from two critical 
problems; fine tuning and coincidence. The first one relates the small value of cosmological 
constant and the second problem caused because we live in an epoch that the magnitude 
of dark energy and dark matter are comparable. Other candidates of dark energy are for 
examples quintessence [1], phantom [2], quintom [3], tachyon [4], holographic dark energy 
[5], K-essence [6] and various models of Chaplygin gas. The simplest case of this model 
based on Chaplygin equation of state [7] to describe the lifting force on a wing of an air 
plane in aerodynamics. The Chaplygin gas (CG) was not consistent with observational data 
[8-11]. Therefore, an extension of CG model proposed [12, 13], which is called generalized 
Chaplygin gas (GCG). However, observational data ruled out such a proposal. Then, GCG 
extends to the modified Chaplygin gas (MCG) [14]. There is still more extensions such as 
generalized cosmic Chaplygin gas (GCCG) [15], and modified cosmic Chaplygin gas [16]. 
On the other hand bulk viscosity plays an important role in the evolution of the universe. 
The idea that Chaplygin gas may has viscosity first proposed by the Ref. [17] and then 
developed by [18-21]. 

In the Ref. [22] a model of varying Chaplyagin gas which interact with a Tachyonic matter 
considered in framework of general relativity, and investigated the potential V{(j)) and ujtot 
numerically. Also, in the Ref. [23] a model of varying generalized Chaplyagin gas which 
interact with a Tachyonic fluid considered and V{(fi) of the Tachyonic fluid and ujtot of the 
mixture investigated graphically. 

In this paper we are going to extend Refs. [22] and [23] to complete version of Chaplygin gas 
and work on the interaction between the mixture of dark matter and dark energy to study 
their behaviors during evolution of the universe. To achieve this purpose sign changeable 
interaction is used to include all the eras from the early time until the end. By numerical 
methods and graphs, the case of varying and the presence of bulk viscosity and cosmic 
elements are considered. Constants are fixed by imposing the constraint V ^ 0, and the 
ojtot implies the quintessence or phantom behaviors. 

We work in the general relativity framework and in FRW metric since has a spatially flat 
homogeneous and isotropic universe that is suitable for Chaplygin gas as a fluid in it and 
cause the acceleration expansion of the universe. 

This paper organized as the following. In next section we review evolution of Chaplygin gas 
(CG) to viscous varying modified cosmic Chaplygin gas (VIVAMCCG) and in section 3 we 
introduce field equations include interaction. In section 4 we study potential V{(j)) of the 
Tachyonic fluid and Utot of the mixture for general case of interaction and various models of 
scale factor. In section 5 we investigate stability of our model and finally in section 6 we 
give conclusion. 



2 C haply gin gas 

In this section we consider perfect fluids as a simple cosmological model and CG [24] is 
the simplest case to fill the expanding universe. The reason to choose CG is based on the 
observational data that the equation of state parameter for dark energy can be less than -1. 
Phantom fields are the same but they have instabilities [25]. The CG equation of state is 
given by the following [26], 

B 

P=--. (1) 

where i? is a positive constant. The CG is also important in holography [4], string theory [27], 
and supersymmetry [28]. Equation (1) is connected to string theory and can be achieved by 
the D-branes Nambu-Goto action which is moving in a (d+2)-dimensional space-time in the 
light-cone parametrization [29] . CG is the only kind of fluid that accepts a supersymmetric 
generalization [30]. Another remarkable characteristics of CG is their Euler equations with 
a large group of symmetry which cause their integrability. CG was studied before [31] in 
the stabilization of branes [32] and black hole bulks [33] . A negative constant B described 
in wiggly strings which causes an anti-Chaplygin state equation [34, 35]. CG describes a 
transition from a decelerated cosmological expansion to the present cosmic acceleration and 
perhaps submit a deformation of A-CDM models. The inhomogeneous CG can combine dark 
energy and dark matter and plays the unification role of them [36, 37]. For studying more 
about CG see Refs. [36-39]. It is also possible to study FRW cosmology of a universe filled 
with generalized Chaplygin gas (GCG) with the following equation of state [40-42], 

B 

with < a < 1. The GCG is also interesting from holography point of view [43]. As we 
can see the GCG is corresponding to almost dust {p = 0) at high density which is not agree 
completely with our universe. Therefore, modified Chaplygin gas (MCG) with the following 
equation of state introduced [44, 45]. 

B 

where /i is a positive constant. This model is more appropriate choice to have constant 
negative pressure at low energy density and high pressure at high energy density. The 
special case of /i = | is the best fitted value to describe evolution of the universe from 
radiation regime to the A-CDM regime. 

The next extension performed by the Ref. [15] where the generalized cosmic Chaplygin gas 
(GCCG) introduced by the following equation of state. 
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This model can also extend to varying modified cosmic Chaplygin (VAMCCG) gas with the 
following equation of state, 
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where uj is corresponding cosmic parameter. Here, B{a) is no longer a constant and depend 
on scale factor via the following relation. 



with, 
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It is also possible to include bulk viscosity. 
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where q is the viscosity coefficient and H = d/a is the Hubble parameter. The equation 
of state (8) corresponds to viscous varying modified cosmic Chaplygin gas (VIVAMCCG) 
which is our interesting case in this paper. 

The total energy density and pressure relate to the mixture of VIVAMCCG interacting with 
Tachyon fiuid and we define, 

Ptot = Ptf + Pde, (9) 

Ptot=PTF+PDE, (10) 

where ptf and ptf denote energy density and pressure of Tachyonic fiuid respectively. Also, 
Pde and pde are dark energy density and pressure which given by the equation of state (8). 



3 Field equations 

As we know the Friedmann-Robertson- Walker (FRW) universe in four- dimensional space- 
time is described by the following metric. 



ds'^ = -dt^ + a{tf 



dr'^ 



kr"^ 



r'dn' 



:n) 



where (if2^ = dO'^ + sin^ 



and a(t) represents the scale factor. The 9 and parameters 



are the usual azimuthal and polar angles of spherical coordinates, with < 6 < tt and 
< (j) < 2n. The coordinates {t,r,9,<p) are called co-moving coordinates. Also, constant k 
denotes the curvature of the space. The curvature k may be not only positive, corresponding 
to real finite radius, but also zero or negative, corresponding to infinite or imaginary radius. 



The possibilities are called closed {k = 1), flat {k = 0), and open (A; = —1). In that case the 

Einstein equation is given by, 

1 
R^lu — T^g^ivR — T^i^, (12) 

where we assumed c = 1, SiiG = 1 and A = 0, k = 0. 

By using the above relations one can obtain the following field equations, 

and 

- = --piptot+Ptot)- (14) 

a o 

The energy-momentum conservation law obtained as the following, 

Ptot + Kptot + Vtot)H = 0, (15) 

Now, for the interaction between dark energy and matter, splits the equation (15) into the 
following equations, 

Pde + ^{Pde+Pde)H = Q, (16) 

where pde is given by the relations (1) to (8) which is Pvivamccg for our complete model, 
and, 

Ptf + Kptf + Ptf)H = -Q, (17) 

where Q is interaction term. Then, the energy density and pressure of a Tachyonic fluid will 
be, 

PTF = 3i/' - PDE, (18) 

and 

-Q - Ptf /-i^n 

Ptf = ^ Ptf, (19) 

where we used the equation (13) and (17). Then, the equation of state parameter and 
potential of Tachyonic matter are given by [22, 23], 

iOTFit) = ^^, (20) 

Ptf 



V{(f)) = V-PtfPtf, (21) 

and Utot is the equation of state parameter of the mixture which given by, 

^„ = ?Il±2SE. (22) 

Ptf + Pde 

In order to investigate potential (21) and the total equation of state parameter (22) we use 
conservation equation (16) to obtain pde- Then, we use the equation (13) to obtain ptf- 
Therefore, the equation (14) gives us Ptf, and hence V{(f)) and utot will be obtained. 



4 Sign changeable interaction and scale factors 

One of the ways to solve the cosmological coincidence problem is to consider the interac- 
tion between the components on phenomenological level. Generally, interaction could be 
considered as a function of energy densities and their derivatives: Q{pi, pi, ...). Interac- 
tions considered as: Q = SHbpm, where fe > is a coupling constant. The general form is 
Q = SHbjpi + 7/9i, where i = TF,DE,tot, and 7 is dimensionless constant [46]. These kind 
of interactions are either positive or negative and cannot change sign. A sign-changeable 
interaction is of the following form, 

Q = qhp + 3bHp), (23) 

where q is the deceleration parameter which is given by, 

^ = -(l + :^)- (24) 

Now, by putting the sign-changeable interaction and the equation of state for viscous varying 
modified cosmic Chaplygin gas (8) in the equation (16), we get to a differential equation as 
the following. 
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H = q{-fp + 3bHp), (25) 



which can be solved numerically for the p during the time. Then p^F and ptf can be 
achieved by equations (18) and (19). Here, the variables are fixed to Vtf -^ 0. 
According to the accelerating expansion of the universe, we work on three kinds of scenarios 
which are based on different eras in the evolutionary process of the universe which all of 
them consist of a kind of expanding exponential scale factor as the foUowings. 

4.1 Emergent scenario 

The universe in the emergent scenario has some interesting properties as the following. The 
universe at large scale is isotropic and homogeneous and there is no time-like singularity. 
The universe is accelerating and may contains exotic matter such as CG. The scale factor in 
this scenario is given by, 

a{t) = ao{B + e^'T, (26) 

where ao > 0, K > 0, i? > 0, and m > 1 [36]. 

First, we consider the simplest case of modified cosmic Chaplygin gas interacting with 
Tachyon fluid in emergent era and reconstruct the tachyon potential and total equation 
of state parameter numerically (Fig. 1). 

By imposing the bulk viscous parameter ? 7^ and variable parameter B{a) in the equation 
of state, we have the VIVAMCCG, which is the general form. Here we suppose various mod- 
els and find their behaviors in respect to their Utot and according to the constraint Vtf -^ 0. 



Analysis of Tachyon potential V{(j)) and Utot shows that in all models for the emergent sce- 
nario and during whole evolution of the universe from the beginning to the end, the total 
equation of state parameter is Utot < — 1 which means the phantom-like behavior (Figs. 
2- 7). This shows that the Emergent scenario is useful for describing the universe because it 
is in agreement with the real data. In the presence of cosmic element in the VIVAMCCG 
in comparison to the VIVAMCG, the Tachyon potential has a faster decreasing in time. In 
VIMCCG and VAMCCG in comparison with the case of MCCG, The presence of viscosity 
and varying condition cause slower vanishing for the potential. 

4.2 Intermediate scenario 

For the scale factor corresponding to intermediate scenario we have, 

a{t) = e''', (27) 

where A > and < /3 < 1 [37, 38]. 

Numerical research on the Tachyon potential and total equation of state parameter in this 

era with the intermediate scale factor reveals the quintessence-like behavior which is based 

on the Utot > —1 (Figs. 8-13). 

In this case, the Tachyon potential of the VIVAMCCG vanishes slower than the VIVAMCG. 

It the VIMCG model, it takes the longest time to Vtf -^ 0. Here, we found a different 

reaction and Vtf in the VIVAMCCG decreases slower by the cosmic effect. 

4.3 Logamediate scenario 

The Logamediate scenario of the universe is motivated by considering a class of possible 
cosmological solutions with indefinite expansion. In this model the scale factor showing the 
accelerating expansion of the universe is given by, 

a(t) = e"(i"W)', (28) 

where a; > 0, and /? > 1 [37]. 

The same process for the Logamediate scale factor gives the Figs. 14-18 which are similar 
to the Intermediate case consistent with the quintessence-like behavior. Here, VIVAMCCG, 
acts has a rather slower decrease in the value of Vtf than the VIVAMCG potential. 

5 Stability 

Now, we are interested to check the stability to find the best model for explaining the 

accelerating expansion of the universe. Positive and bounded squared sound velocity of 

Chaplygin gas is a remarkable feature that is a non-trivial fact of fluids based on their 

negative pressure 

dP 



Here, we study the stability of these variety of models with respect to the positive quantity 
for v"^. This is important to mention that the parameters are fixed according to Vtf — ;■ in 
the interaction. 

5.1 Emergent scenario 

In this case, the models VAMCCG, VAMCG, VIMCG and VIVAMCG are stable during the 
whole time and the model of MCCG has no stability here. For the model of VIMCCG, there 
is no stability until t = 15 and for VIVAMCCG, as the general form, the stability is present 
from the beginning but after the time t = 108, no stability exits in the model. 

5.2 Intermediate scenario 

In the general model of VIVAMCCG the stability stays only till early moments t = 9 and 
stability continues no more. That is like the emergent scenario but with this difference 
here that the stability vanishes so soon in comparison with the similar model in emergent 
kind. The models VIMCG, VAMCG, VIMCCG and VIVAMCG are stable all the time of 
evolution. In the VAMCCG the stability disappears very soon at t = 3 that is a distinct 
situation of this scenario. 

5.3 Logamediate scenario 

For VIVAMCCG in logamediate scenario, the stability holds in early time before t = 18 and 
exactly like the other scenarios during the late time the model misses stability. In the other 
models as VIVAMCG, VAMCCG and VAMCG the stability is always available during the 
whole evolution of the universe in respect to the time. But in the VIMCCG as the emergent 
similar case, the stability holds until the time t = 432. 

6 Conclusion 

In this paper we extend Ref. [22] and [23] to the case of viscous varying modified cosmic 
Chaplygin gas and studied interaction with Tachyonic fluid. We considered general form of 
interaction and three different cases of the universe which are Emergent, Intermediate and 
Logamediate scenarios. In the Emergent case, we found the phantom-like behavior due to 
the ujtot and quintessence-like behaviors in the Intermediate and Logamediate era is obvious. 
The comparison between the values of cjtot in the case of MCCG and in the presence of 
viscosity shows a tiny difference that is negligible. Also, The cosmic element causes no 
comparable difference in the Utot- By taking <;" = 0.3 the values of Vtf and Utot decreases a 
little. 

The viscosity causes a little delay in the vanishing procedure of the Tachyon potential and 
the cosmic element appears by the same feature. 

Stability of the general model of VIVAMCCG of all the three scenarios are contemporary 
during the time and vanishes in the late time of expanding evolution of the universe. It 
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shows that the viscous varying modified cosmic Chaplygin gas model is vahd for the early 

universe. 

Since the Phantom-like behavior is more consistent with the observational data, Then the 

Emergent scenario is better than the Logamediate and Intermediate ones and the stable 

models in this scenario are the best ones among all the cases we have considered here. These 

suitable models are VAMCCG, VAMCG, VIMCG and VIVAMCG in Emergent scenario. 
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Figure 1: Plot of V and Utot in terms of time with m = 1.1, A 
B = 1, ao = 0.6, n = 0.3, u = -0.5, a = 0.5 and K = 0.03. 
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Tachyon potential 




(Emergent): VIMCG 

Figure 2: Plot of V and Utot in terms of time with m 
ao = 1, yu = 0.3, ^ = 1, a = 0.5 and K = 0.02. 
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Figure 3: Plot of V and Utot in terms of time with m = 1.1, b = 0.1, 7 = 0.7, B = 1, ao = 0.6, 
i = -0.6, j = -0.08, /i = 0.3, a = 0.5 and K = 0.03. 
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Figure 4: Plot of V and Utot in terms of time with m = 1.1, b = 1, K = 0.03, 7 = 0.7, B 
ao = 0.6, i = -0.6, j = -0.8, /i = 0.3, a; = -0.5 a = 0.5 and K = 0.03. 
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Tachyon potential 



Total EOS parameter 




(Emergent): VIMCCG 



(Emergent): VIMCCG 



Figure 5: Plot of V and Utot in terms of time with m = 1.1, A = 0.1, b 
7 = 0.7, B = l,ao = 0.6, /i = 0.3, ^ = 1, a; = -0.5, a = 0.5 and K = 0.03. 
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Figure 6: Plot of V and Utot in terms of time with m = 1.1, b = 0.1, K = 0.03, 
S = 1, ao = 0.6, i = -0.6, j = -0.08, iJ = 0.3, ^ = 1, a = 0.5 and K = 0.03. 
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Figure 7: Plot of V and Utot in terms of time with m = 1.1, b = 1, K = 0.03, 7 = 0.7, B = 1, 
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Figure 8: Plot of V and Utot in terms of time with (j, = 0.3, 7 = 0.7, A = 0.7, /3 = 0.8, >? = 1, 
A 



1,6 



0.1, and a = 0.5. 
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Figure 9: Plot of V and ujtot in terms of time with /i = 0.3, 7 = 0.7, A = 0.7, j3 = 0.8, 
b = 0.1, i = -0.6, j = -0.03 and a = 0.5. 
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Figure 10: Plot of V and Utot in terms of time with /i = 0.3, 7 = 0.7, A = 0.7, /3 = 0.8, >? = 1, 

-0.5 and a = 0.5 



A = OA, 6 = 0.1, u 



16 



Tachyon potential 




CO 





Total EOS parameter 


0.88- 




0.90- 




0.92- 




0.94- 




0.96- 


\ 


0.98- 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



(Intermediate): VAMCCG 



10 20 30 40 50 60 70 80 90 100 

t 

(Intermediate): VAMCCG 



Figure 11: Plot of V and Utot in terms of time with /i = 
i = -0.6, j = -0.03, b = 0.1, CO = -0.5, and a = 0.5. 
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Figure 13: Plot of V and Utot in terms of time with ji = 0.3, 7 = 0.7, A = 0.7, /3 = 0.8, >? = 1, 
i = -0.6, j = -0.03, b = 0.1, u = -0.5, and a = 0.5. 
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Figure 14: Plot of V and Utot in terms of time with x = 0.05, b = 0.4, /3 = 2, 7 = 0.1, 
a = 0.5, /i = 0.3, i = -0.8 and j = -0.4. 
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Figure 15: Plot of V and Utot in terms of time with x = 0.1, b = 0.4, /3 = 2, 7 = 0.1, a = 0.5, 
1^ = 0.3, u = -0.5, i = -0.8 andj = -0.4. 
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Figure 16: Plot of V and Utot in terms of time with x 
a = 0.5, fj, = 0.3, <;■ = l,u = —0.5 and A = 1. 



0.05, b = 0.4, /3 = 2, 7 = 0.1, 
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Figure 17: Plot of V and Utot in terms of time with x = 0.5, b = 0.4, /3 = 2, 7 = 0.1, a = 0.5, 
^ = 0.3, q = l,i = -0.8 and j = -0.7. 
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Figure 18: Plot of V and Utot in terms of time with x = 0.5, b = 0.4, /3 = 2, 7 = 0.1, a = 0.5, 
/i = 0.3, ^ = 1, z = -0.8, j = -0.7 and u = -0.5. 
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